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Abstract 


We consider the maximum flow network interdiction problem. We provide 
a new interpretation of the problem and define a concept called ” optimal 
cut”. We propose a heuristic algorithm to obtain an approximated cut, and 
we also obtain its error bound. Finally, we show that our heuristic is an 
q@-approximation algorithm for a class of networks. By implementing it on 
three network types, we show the advantage of it over solving the model by 
CPLEX. 
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1 Introduction 


In the maximum flow network interdiction problem (MFNIP), an interdic- 
tor selects a subset of arcs in a given capacitated network to fully interdict 
them subject to a given budget limitation in order to minimize the evader’s 
maximum flow in the remaining network. Two positive numbers are assigned 
to all arcs as their capacities and interdiction costs. The problem is known 
as cardinality maximum flow network interdiction problem (CMFNIP), if all 
interdiction costs are equal to one. 

The original form of MFNIP roots back to more than five decades ago, 
(see [25, 10]), and has followed by researchers with several variations and 
application fields. Wood [26] provided the first 0-1 mathematical program- 
ming with two types of valid inequalities for MFNIP and showed that the 
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problem is NP-hard (Nondeterministic polynimial time), even in cardinality 
case. He also added two types of valid inequalities to the relaxation of the 
problem both in general and cardinality cases. Altner, Ergun, and Uhan 
[3] suggested two types of polynomially separable inequalities for CMFNIP 
and showed that even if the known IP model of the problem defined in [26] 
is strengthened, with these valid inequalities, there is still large integrality 
gap. Two extended formulations are provided for CMFNIP defined in [2], 
one based on the valid inequalities of [3], and the other based on some new 
valid inequalities. 

In a stochastic interdiction problem, some or all parameters are consid- 
ered as stochastic variables (see [9, 19, 24]). The maximum reliability path 
problem on bipartite networks in the stochastic case is the subject of [23]. 
They provided a convex hull for the polytope of interdictors decision vari- 
ables and strengthened the relaxation model up to 25% by adding separable 
inequalities. If there exist information asymmetries between interdictor and 
evader, then a more computationally difficult version of the problem will arise 
(see [14, 21]). 

Interdicting arcs on dynamic network, in which a positive number is also 
assigned to all arcs as their traversal time, was studied in [1]. They minimized 
the dynamic maximum flow in a given period of time. There is also another 
dynamic version in [13], where interdiction strategies may vary in discrete 
periods of time. 

Network interdiction problems have several applications such as the in- 
fection control in hospital (see [4, 17]), identifying critical infrastructure (see 
[8, 5]), nuclear smuggling (see [16]), vulnerability analysis (see [5]), supply 
chain networks (see [12]), and border control (see [19, 22]). Minimizing the 
connectivity of a given network by interdicting a path is another NP-hard 
interdiction problem, called the critical disruption path; see [11]. 

There exists only a (1+e,1+ +)-pseudo polynomial algorithm for MFNIP; 
see [6]. For planar network, there is a fully polynomial time approximation 
scheme (FPATAS) for it; see [20]. Recently Chestnut and Zenklusen [7] pro- 
vided a 2(n — 1)-approximation algorithm for general graphs with n vertices. 

In this paper, we apply the LP (Linaer programming) relaxation of the 
model defined in [26] and provide a heuristic algorithm for MFNIP. We also 
calculate the error bound of the method in gencral networks, while for some 
special networks with proven large integrality gap, we will show that the 
proposed algorithm is an a-approximation method. 

This paper is organized as follows. In Section 2, we consider the math- 
ematical model for MFNIP and provide a new interpretation for MFNIP. 
Also we define the concept of optimal cut”. An st-cut is recommended in 
Section 3 as an approximation for the optimal cut, which leads to a heuristic 
method. The error bound of the heuristic algorithm is calculated in Section 
4. Moreover, we show that the proposed heuristic is an a-approximation al- 
gorithm for a class of networks, which overcomes the algorithm defined in [7] 
for this special class. Section 5 reports the numerical results for our approach 
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on three types of networks, and we compare the performance with the 0 — 1 
model in literature. Finally, Section 6 concludes the paper. 


2 Interpreting MFNIP 


Let G = (N,A) be a digraph with set of nodes N and set of arcs A. Let 
N = (G,u,r, s,t) be a given network, in which u and r are the vectors of arc 
capacities and interdiction costs, respectively, and {s, t} is the set of terminal 
nodes. Our goal in MFNIP is to find the subset of arcs with interdiction costs 
lower than or equal to a given limited budget R, whose removal leads to the 
least maximum st-flow in the remaining network. 

The first mathematical model for MFNIP was proposed in [26], with the 
following decision variables: 


1 if node 7 is in the sink side of the cut, 
T= 
: 0 otherwise, 


9 f if arc (1,7) is the forward arc of the cut, 
ig = 


O otherwise, 


8 f if arc (i,j) is the forward arc of the cut and not interdicted, 
y= 


0 otherwise, 


if arc (2,7) is interdicted, 
Vig = 


0 otherwise. 


Consider the following integer programming from [26]: 


min oy Uig Biz, (1) 


(4,7)E€A 
st. 1 —1j +6; >0, (1,9) € A, (1a) 
Tm —Ts, > 1, (1b) 
:S Tig Vig <A, (1c) 
(i,5)EA 
Vet he-lg 20 Wj) eA, (1d 
m€E{0,1}, tEN, (le 
Oi; € {0,1}, (7) € A, ( 


Vig © 20, A (2,9) € A, ( 
Big € qd fs (7,7) eA, (1h 
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In optimality, inequality (1d) changes to the equality (see [26]). Hence, 
letting 6;; = 6i; + yi; leads to the main integer programming defined in [26]. 
Here we substitute 8;; by 0:; — yj to obtain 


min )) tay(Bi3 — 4s) (2) 
(4,j)EA 

at. Yay < O33, (2,9) € A, (2a) 
(1a) — (le), 
m €{0,1},7E N, 
Vig, 917 € {0,1}, (2,9) € A. 


Note that 6;; € {0,1}, for all (,7) € A, so inequalities (2a) are imposed 
to model (2), since no arc out of the cut will be interdicted in optimal case. 
The objective function (2) can be rewritten as follows: 


min ; Ujj9;; — max s Wig Vig 
T, . y : 
(i,j)EA (t,j)€A 


which means we are solving the minimum capacity cut problem with con- 
straints (la), (1b), (le), and (1f) and the knapsack problem with constraints 
(1c) and (1g), and constraints (2a) relate these two problems. 


Remark 1. Solving the 0 — 1 knapsack problem on a cut, means to sup- 
pose forward arcs of the cut as items, arc capacity as the item’s value, and 
arc interdiction cost as the item’s weight. The limited budget R is also the 
limited capacity of the knapsack. Interdicting arc (7,7) corresponds to select 
item (i,7) for knapsack. 


Remark 2. The MFNIP is to solve the 0 — 1 knapsack problem on all 
feasible st-cuts of the problem and choose the cut, in which the remaining 
capacity is minimum after interdiction. 


Therefore we define the following concept. 


Definition 1. For a given network NV = (G,u,r,s,t¢) and budget R, an 
st-cut C is called optimal cut, if solving 0 — 1 knapsack problem on it (based 
on Remark 1), results in minimum remaining capacity on it. The optimal 
cut is shown by C*. 


Now it is obvious that MFNIP is to find the optimal cut according to 
Definition 1, which is not necessarily the minimum capacity cut. 
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3 Approximated cut 


In order to approximate the optimal cut, consider the dual of the LP relax- 
ation of model (1) as follows: 


max f— Rw (3) 

at SO fem YS eae (3a) 
I:(8,j)EA T:G,8)EA 

> Jee SHH (3b) 
F:4G,9)EA PGDEA 

», feo SD) fe=-, (3c) 
DL IVEA TG t)EA 

Aag = az, (G9) © Ay (3d) 

Aap — wWreg <0, (4,9) € A, (3e) 

fg se sl. Ga) eA; (3f) 

Vij fig 2 0, (4,9) € A, (3g) 

f,w > 0. (3h) 


The variables f;; and f are the dual variables of constraints (1a) and 
(1b), respectively, and the variables \;; and w are also the dual variables of 
constraints (1d) and (1c), respectively. 

For a fixed w, model (3) is a maximum flow problem reduced by a factor 
Rw, in which the capacity of arc (7,7) is 


Vij = min{u,;, wry}. (4) 


Let w* be the optimal value of w in model (3). Then for w = w%*, the 
maximum flow value is equal to f* — Rw* and is obtained by capacities 
Aig = min{u,;,w*ri;}. The dual problem of maximum flow problem (3) is 
the following minimum capacity cut problem: 


Uc,,. = min » rig Oi (5) 
(4,j)EA 
Tt, — 7; +6;; >0, (4,9)E A (5a) 
nan St (5b) 
6:5 > 0, (4,7) € A. (5c) 


Arc capacities in model (5) are calculated by (4) with w = w*. Also C,,+ 
is the cut obtained by model (5) and Uc,,,. is its capacity. 
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Remark 3. Since w is the dual variable of the LP relaxation problem, 
the equality C,,» = C* does not necessarily hold. 


3.1 Parametric procedure 


In this section, we provide a parametric procedure to calculate w*. For a 
fixed w, let f(w) be the optimal value of the maximum flow problem in (3), 
according to arc capacities in (4). It is clear that f(0) = 0 and that increasing 
w increases f(w) up to w = Wmax = MaX(i,j;,e4 =, SO 
ig 
max f(w) = FW) = fmax- 
Let E(w) be the set of arcs for which their capacities are decreased by the 
parameter w: 
E(w) = {(4, 9) € Alig = wrig}- 

In this paper, E(w) is called set of active arcs. 

If Uc shows the capacity of s — t cut C, f(w) is the maximum flow, and 


Z(w) is the optimal objective value of model (3) by arc capacities obtained 
by (4), then by the maximum flow-minimum cut theorem, we have 


f(w) = Ue, = > Wij + Ww 3 rij. 


(4,7) €Cw\ E(w) (4,7) €Cw E(w) 


It is not difficult to verify that f(w) is a concave piecewise linear function, 
since increasing w leads to decrease number of arcs in E(w) and consequently 
in Cy, E(w), so increasing w decreases the slope of linear pieces. Let 


Z(w) = f(w) — Rw = 1 Uy + w ( es rg = Hy (6) 


(i.9) €Cw \ E(w) (4,9) €CuN E(w) 


Then for w € [0, fsx), Z(w) is a concave piecewise linear function as 
shown in Figure 1. The slope of any linear piece of Z(w) is 


s(w) = S- rig — R. (7) 


(i,j) €CwN E(w) 


Note that similar to f(w), increasing w limits the set E(w) and the slope 
s(w) decreases, so Z(w) is a concave function. 

Now in order to calculate Z(w), it suffices to start with any w; and w, 
with 0 < w; < w, and obtain the cross point of the line passing through point 
(w;, Z(wz)) with slope s(w;) and the line passing through point (w,, Z(w,)) 
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with slope s(w,). Call this point (w,Z). Whenever Z = Z(w), we stop, 
otherwise we update w; and w,. Algorithm 1 formally explains the procedure. 


Algorithm 1 Parametric procedure 
Input: G = (N, A), NV = (G,u,r, s,t); Total budget: R. 
Output: Optimal value w*, for problem (3). 
1: Let we = 0 and w, = maxi, jc ae 
2: For w = w; and w = w,, calculate Z(w) and s(w) from (6) and (7), 
respectively; 
3: Let wy = 2d Zlue we sure) we se) and Z = Z(we) + (w — we)s(we); 
4: Calculate Z(w) and s(w); 
5: if Z(w) = Z then 
6 w=; 
7 
8 
9 


if s(w*) > 0 then 
Let (wz, we) = (w*, wr); 


else 
10: Let (w7, wz) = (we, w*); 
11: end if 
12: Stop; 
13: else 
14. if s(w) > 0 then 
15: Let w, = ©, Z(we) = Z(w) and s(we) = s(w) ; 
16: end if 
17; if s(w) <0 then 
18: Let w, = ©, Z(w,) = Z(w) and s(w,) = s(w) ; 
19: end if 
20: end if 


21: Go to step 3 


In step (1), lower and upper bounds of w are determined. In step (2), 
values of Z(w) and s(w) are obtained. In step (3), the crosspoint of two 
lines passing from (we, Z(we)) and (w,, Z(w,)) with slopes s(wz) and s(w,), 
respectively, is achieved. In step (4), values of Z(#) and s(w#) are obtained 
by using equations (6) and (7), respectively. In step (5), we check the stop 
condition, if it is not satisfied, and we update values of wz and w,, in step 
(13) and repeat the procedure. 


Note that calculating Z(w) and s(w) for a given w in steps (2) and (4), 
means to solve a maximum flow problem for which there exist several al- 
gorithms, for example, the algorithm of [18], which has a running time of 
o(|N||A|). Therefore, Algorithm 1 is running in polynomial time. 
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Table 1: Notations 


ZIP Optimal objective value of model (1). 

Zrip Optimal objective value of model (3). 

Zu Objective value obtained by applying Algorithm 2. 
f(w) Maximum flow with capacities (4). 

Z(w) = f(w) — RW. 


Eo = {(4, 9) €E ie Ui; < w* rij}. 

Ey = {(4,J) € Cyx|uij > w* ris}. 

Ty = {(i,j) € ae |(t, 7) is interdicted by Algorithm 2, 
(or is selected in knapsack problem) }. 

Io = {(4, 3) € Iu ij Swag}. 

qi = {(t,7) € Iu+|uiy > w* ri, }- 


3.2 Heuristic method 


Now we suggest a heuristic method, and in section 4, we will show this heuris- 
tic turns out to be an a-approximation method for a class of networks. 

The proposed heuristic is based on the cut C,+«. In this method, we con- 
sider cut C,,« as an approximation for the optimal cut C*, and we call it an 
approximated cut, then we solve knapsack problem on it. 


Algorithm 2 The heuristic algorithm 
Input: NV = (G,u,r,s,t), G = (N,A); Total budget: R. 
Output: Subset of arcs of cost lower than or equal to R 

1: Calculate w* by algorithm 1; 

2: Obtain st-cut C,,«, from model (5); 

3: Solve 0 — 1 knapsack problem on C,,«, according to Remark 2. 


It is clear that Algorithm 2 is not polynomial because of step (3), how- 
ever one can apply some known heuristics, for example, greedy algorithm to 
approximately solve the respective knapsack problem. 


4 Error term 


In this section, we calculate the error term of Algorithm 2. Table 1 shows 
notations, applied in this section. 

First, note that f(w) is a linear function in of w and Z(w) is a concave 
piece-wise linear function of w, starting from original and meets the vertical 
axis in w = ees where fmax is the maximum value of f(w) among all values 
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of w. For a fixed w, the slow of any piece of Z(w) is calculated as follows: 
s(w) = a tig — RK. 
(J) EC*N{G,F) rig <uay} 


Figure 1 shows the function Z(w), where wy and w,. are two arbitrary points 
lower and greater than w*, respectively. 


Z(w) 


Figure 1: Diagram of Z(w) 


According to Table 1, we have C,,» = Eo U Ey and I,» = Ip UL,. It is 
clear that 


Zrip = f* —w*R=U¢,. —-w R= S> Uiy + w*( » riz — R). (8) 
(4,9) € Eo (¢,j)€E1 


Also 


ZH = ‘S Uizg — > Uig 


(4,5) €Cw* (4,5) €Iw* 
S Do ous DO uy (9) 
(4,9) €Cw* (,9)€h 
< > Wij — Ww" > liz: (10) 
(4,5) €Cw* @7)Eh 


Now it is easy to conclude 
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Zrip <Zrp <Zy 
=> Zy —-Zip <Zy — Zrip 


< y Wig — Ww ) rij — ) Uij — W™ S rig tw'R 


(I) ECw* (jen te Eo (i j)EFA 
= S> Uiy — Ww" ys rig tw (R- >> Tig): (11) 
(45) EF 1 (EN (4,5) EF 


Therefore the relative error is 


Zy — Zp Z Zy — Zip 


Err = < (12) 
ZIP ZRIP 
Z ER wig — w* wen rig tw'(R- GER rij) 
7 i,j) eB Vis + W* (Di. jyen, Tis — FR) 
LG DEF, wig -w" LG jer . x 
4 ; ifs(w*) > 0, 
(i,j) € Eg Mis f ( ) 
a (13) 
Da wig wen rig tw"(R-Via nem ri) if s(w*) <0. 


Die By Vii FW* QUE New, THF) 


4.1 Approximation scheme for I,,,, networks 


The J,,,,, class of networks (see [3]) is defined as follows: 

For two positive integers p and &, with k > 2 and wp >> k, Inj = (N, A), 
where N = {s,t} UX UYUZ, A= ASUAPUALUA?, with |X| = «, 
IY] =|Z|=y, and 


A® = {(s,%), 7€ X}, Uz =e, (4,7) € A’, 
At — {(2,¢), LEX UY}, Uj = 1, (i, 7) E At, 
A’ ={(i,j), 1€ XG EY}, uiyg=p?, (i,7) € A®, 
At = {(s, i), (i,t), 1 € Z}, wy =u, (i,9) € At. 


All interdiction costs are equal to one and R= w+«-—1. Altner, Ergun, 
and Uhan [3] improved valid inequalities of [26] and showed that adding their 
valid inequalities does not improve the value of the objective function. They 
also proved that for J,,,, networks, the integrality gap of model (1) is not 
bounded by a constant. According to Definition 1, {(s,i), i © X UZ} is 
the optimal cut (see [3]). Solving the knapsack problem on the optimal cut 
results in Zy7p = yw amount of remaining capacity. 

In this section, we show that Algorithm 2 finds a bound for integrality 
gap of J,,,,, networks. 
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Figure 2: Network instance of Ix, with capacities 4;; = min{u;;,w*} 
and p< w* < 1? 


Lemma 1. Let w* be optimal value for model (3) for an instance of J, ,,; 
then w* < p. 


Proof. Note that 


w* <wWmax = max HD iP 

(I)EA Vij 
Now by contradiction suppose that ps < w* < px2; then we have the network 
in Figure 2. All arcs in Figure 2 are labeled with two numbers (a, |b|), where 
a is the arc capacity according to (4), and b is the number of such arcs in the 
network. It is clear from Figure 2 that {(i,t), 7€ X UY UZ} is the minimum 
capacity cut with capacity « + u+w*p. Consequently the optimal value of 
the objective function of model (3) is equal to 


Z(w*) =pt+K+w*p— Ru. 


Comparing two cuts with arc capacities in (4), one with w = w* and one 
with w = p, leads to 


Z(u) — Z(w*) = pw +p +e— Ru-p—K—w*p+ Ru* 
=p eh RO =p) 
= —w*pt (ut+K—1)(w* — p) 
=kKw* —w* — pkK+p=(w* —p)(K-1) >0 
= Zw") < 2(u), 


which contradicts with optimality of w*. Therefore w* < wp. 


Theorem 1. Algorithm 2 on an instance of J,,,,, is 1_approximation algo- 
, le 
rithm. 


Proof. Consider an instance J,,,,, and assume that w* has been obtained by 
model (3). Now we have to solve the minimum capacity cut problem on 
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Figure 3: Network instance of I,,,, with capacities A;; = min{u;,;, w*} 


the network with arc capacities \;; = min{u,;, w*} for all (i,7) € A. This 
network is shown in Figure 3. 

Note that in arc capacities of Figure 3, we have applied Lemma 1, that is, 
min{p,w*} = w*. First we have to find the minimum capacity cut. Among 
all st-cuts in Figure 3, there are two st-cuts with the minimum capacity de- 
pending on the value of w*: The capacity of dotted cut is equal to hw* +w”* py, 
and the capacity of the solid cut is equal to w+ + w*yu. Therefore we have 
the following two cases: 


Case 1: If «(w* — 1) < yp, then the dotted cut is the min cut obtained 
by model 5, since we mentioned that this cut is the optimal cut and solving 
knapsack on it will give us the optimal value equal to u. 


Case 2: If K(w* — 1) > w, then the solid cut is the result and knap- 
sack problem should be solved on it. We have to remove «+ y — 1 arcs to 
minimize the remaining capacity. By removing all arcs from set Z to node 
t and all arcs but one from set X to t, the remaining capacity is equal to +1. 


Note that if «(w* — 1) = yp either first case or second case may happen. 
Thus 


4H ~ Zip hk +1— pb _ 1 


ZIP Hl be 


5 Numerical results 


In this section, we solve model (1) for three types of networks and also apply 
Algorithm 2 on them to investigate the efficiency of the proposed algorithm. 
We carry out all computational tests on an Intel(R) Core(TM)2 Duo Proces- 
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Table 2: Results for Ic, networks 


Network parameters Optimal value CPU time (min:sec) 
Ll K (\N|, |Al) ZRIP ZH model (1) Algorithm 2 
10 2 (23,54) 6.00 11 00:00.509 00:00.692 
20 5 (46,170) 5.00 20 00:00.405 00:00.666 
40 5 (86,330) 9.00 40 00:00.494 00:00.739 
50 5 (106,410) 11.00 51 00:00.415 00:00.818 
100 =10 (211,1320) 11.00 100 00:00.991 00:01.241 
150 = =20 (321,3490) 8.50 150 00:04.707 00:01.541 
150 =50 (360,8050) 4.00 150 00:14.781 00:02.303 
200 850 (451,10700) 5.00 200 00:24.617 00:03.083 
200 70 (471,14740) 3.86 200 00:45.182 00:04.316 
200 100 (501,20800) 3.00 200 01:25.428 00:06.858 
500 100 (1101,51700) 6.00 500 12:59.496 00:21.992 


sor 2.20 GHz, 4 GB of RAM, GAMS 24.2.2 generates models, and CPLEX 
12.6.0.0. solves them. 


5.1 I,,,,, networks 


I, networks were introduced in section 4.1. Table 2 shows results for this 
class of networks. Eleven J,,,,, networks with different parameter values were 
tested. The column “Optimal value” shows the value of objective function. 
Also Zprp is the optimal value of relaxation of model (1), which shows the 
large integrality gap, while Zy is the value obtained by Algorithm 2. Note 
that there is no column for optimal value Z;p, since we know that Zp = p, 
so this value is not reported in this column. 

It is clear in Table 2 that specially for large amount of « and yp, Algorithm 
2 is extremely faster than solving the original model (1) by CPLEX. However 
results in section 4.1 also confirms the advantage of Algorithm 2. 


5.2 Gp,g networks 


In this section, we test the proposed algorithm on some layered networks 
called Gp,,, with g columns of nodes and h nodes in each column. Nodes of 
the first column are considered as source nodes and nodes of the last are sink 
nodes. Any node in the 7th column is incident to all nodes in the 7 + 1th 
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Ss > ot 
Figure 4: Network instance of Gp, 
column, fori =1,...,g—1. Arc capacities and interdiction costs are random 


integer numbers, uniformly distributed in [10, 30]. Figure 4 shows a schematic 
of Gp,, networks. 


Table 3 shows the results obtained by solving model (1) and applying 
Algorithm 2 on Gg networks. Despite [,,,, networks, Algorithm 2 does not 
perform well on G;,,, networks in compare to model (1). Moreover there are 
several cases in Table 3, where Algorithm 2 does not find the optimal value. 


5.3 Star mesh networks 


This section describes tests on star mesh networks; see [15], for a detailed 
description of this network type. Figure 5 shows a star mesh network. Also 
Sp» is a star network with p rays and v rings. The node in the center 
of the star is the source node s, and all nodes in the outer ray are sink 
nodes. Arc capacities and interdiction costs are drawn form a discrete uniform 
distribution on [1, 10] and [10, 30], respectively. 
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Network parameters Optimal value CPU time —_(min:sec) 

(h, g) (\V|, |A]) R Zrip Zip Zy model (1) Algorithm 2 
(5,5) 27,110 43 99 103 105 00:01.974 00:01.741 
(5,10) 52,235 21 207.2 209 209 00:00.550 00:01.721 
(6,10) 62,336 235 230.23 242 253 00:01.250 00:01.140 
(8,10) 81,592 235 670 673 680 00:03.160 00:01.156 
(7,9) 65,406 75 439.8 448 469 00:01.043 00:01.085 
(9,9) 83,666 82 455.105 461 482 00:01.121 00:01.083 
(9,15) (137,1152) = 251 228.28 239 239 00:02.930 00:01.488 
(10,8) 82,720 257 571.14 574 577 00:00.701 00:01.103 
(10,15) (152,1420) 300 290.60 293 293 00:01.144 00:01.582 
(25,30) (751,8175) 700 5750.75 5751 5760 01:00.500 02:15.684 
(25,30) (751,8175) 2000 1581.25 1582 1585 03:02.430 03:36.746 


Figure 5: Star mesh network instance of Sp,v 


Table 4 provides detailed results on star mesh networks. For fairly small 
networks, model (1) is much more faster than Algorithm 2, however as net- 
work parameters increase, Algorithm 2 exceeds model (1). 


6 Conclusion 


In this paper, we applied the 0 — 1 mathematical model for MF NIP in liter- 
ature and provided an interpretation and consequently a heuristic algorithm 
for it. We also showed that the provided algorithm is a “approximation al- 
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Table 4: Results for star mesh networks 


Network parameters Optimal value CPU time —_(min:sec) 
(p,v) |V|, | Al) R Zrip Zip Ly model (1) Algorithm 2 
(5,50) 251,505) 30 3.42 4 4 00:00.609 00:00.942 
10,10 101,210) 25 15.42 17 18 00:00.339 00:00.932 
15,10 151,315) 30 29.3 31 32 00:00.397 00:00.840 
20,15 301,620) 71 34.27 35 35 00:00.490 00:00.990 
20,20 401,820) 100 29.66 31 31 00:01.176 00:01.033 
25,20 (501,1025) 130 30.96 32 33 00:00.881 00:01.227 
30,30 (901,1830) 200 29 30031 00:02.882 00:01.947 
40,30 (1201,2440) 200 49.3 50 ~—-50 00:04.719 00:03.209 
40,50 (2001,4040) 300 33.54 34 36 00:20.667 00:09.278 
70,50 (3501,7070) 500 54.85 55 55 01:36.782 00:36.621 
70,80) (5601,11270) 1000 6.26 7 7 06:30.792 02:14.322 


gorithm for J,,,,, class of networks. However the error of the proposed method 
was also calculated for general networks. 

Numerical results showed that the proposed algorithm is generally faster 
than the 0 — 1 model by CPLEX, specially on large size grids and star mesh 
networks. Although the proposed algorithm is not polynomial due to the 
knapsack problem in it, but it is still faster than 0—1 model on large networks. 
The reason is that according to the provided interpretation of the problem, 
solving MFNIP is equivalent to solving knapsack problem as many as st-cuts 
in the network, while our proposed heuristic algorithm reduced it to solving 
knapsack problem just once. This fact indeed, explains its preponderant 
performance. 
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